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The B3-VLA CSS sample. II: VLBA images at 18 cm 
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Abstract. VLBA observations at 18 cm are presented for 28 Compact Steep-spectrum radio Sources (CSSs) from 
the B3-VLA CSS sample. These sources were unresolved in previous VLA observations at high frequencies or 
their brightness distribution was dominated by an unresolved steep spectrum component. More than half of them 
also showed a low frequency turnover in their radio spectrum. The VLBA images display in most cases a compact 
symmetric structure. Only in a minority of cases complex morphologies are present. 
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1. Introduction 

The population of CSSs (Compact Steep-spectrum 
Sources) & GPSs (GHz Peaked-spectrum Sources) (see 



O'Dea, 1998 for a review) has received recently an increas- 
ing attention as at least the two sided (or symmetric) ob- 
jects of this class are considered to be the young stage in 



al 



1996), from the Compact (CSO, size 

i 



the radio source evolution (Fanti et al. 1995; Readhead et 

kpc)]J 



< 1 h-> 
< size < 20 h~ 



through the Medium (MSO, 1 h~ L kpc 
kpc) to the Large Symmetric Object (LSO, size > 20 hr 1 
kpc) . A strong support for this interpretation comes from 
the large expansion velocities recently measured i n a h and- 
ful of sub-kpc sources (see compilation in Fanti 200C and 
references therein) and from the determination of the ra- 
diative ages of a larger number of CSSs spanning a range 
of lin ear size from sub-kpc to ~ 20 ft, _1 kpc (Murgia et al. 
1999b . 

by 



1982 



Starting from the seminal paper by Baldwin 
models for radio source evolution have been presented by 
several authors which link the CSS/GPS pop ulation to 
the larger size so urces, (see e.g. Begelman 1996 ; Kaiser & 
Alexander |l997|) . 

All models essentially predict a decrease of the source 
luminosity with time, and therefore with increasing source 
linear size. A key test for the models is the ability to re- 
produce the distribution of sources in the "Radio Power - 
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Linear Size" (Pr-LS) plane (Baldwin |l982j). CSSs/GPSs 
seem to be in the right proportion as compared with model 
expectations (Fanti et al. 1995; Readhead et al. 1996[ ). 



However, from the samples studied so far the statistics 
have been not very large and suspicions were raised on 
some "irregularities" in the LS distribution (O'Dea & 
Baum, 1997) and on the, perhaps too high, number of very 
small size C SSs, the GPSs (Polatidis et al. |1999[ Fanti & 
Fanti |2001[ ). 

Recently we have selected a new sample of CSSs from 
the B3-VLA catalogue (Vigotti et al. |1989| ) aimed at in- 
creasing significantly the statistics on sources with LS in 
the range 0.4 frr 1 < £S(kpc)< 20 /i" 1 (Fanti et al. [2001 



hereafter referred to as Paper I). This sample, consisting 
of 87 CSSs (LS < 20 hr 1 kpc), has VLA observations at 
1.4 GHz (A and C configurations), 5 and 8.4 GHz (both 
A configuration). About 50 % of these sources were not 
resolved or were poorly resolved even at 8.4 GHz (beam 
size ~ 0.2"). For them two VLBI observing projects were 
undertaken with the VLBA and EVN & MERLIN. 

This paper reports on the results of the VLBA observa- 
tions of the most compact sources. The EVN & MERLIN 
observations instead arc presented in a companion paper 
by Dallacasa et al. fl2002|, Paper III). 



2. The Source list and VLBA observations 

Among the twenty-eight sources selected for VLBA ob- 
servations, twenty-four are unresolved or slightly resolved 
(angular size < 0.2") in the VLA 8.4 GHz observations. 
The remaining four sources are dominated in flux density 
by an unresolved or slightly resolved component, with the 
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Table 1. The VLBA Sample: See text for a description of the columns. 



Source 


Id 




z 


LAS 


Log Po.4GHz 


a VLA 
S>1.67 


o VLBA 

& 1.67 


LASvhBA 


LLS 


Morph. 










mas 


W/Hz h~ 2 


mjy 


mjy 


mas 


kpc n 






(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


0039+373 


G 




1.006 


100 


27.37 


784 


779 


120 


0.50 


cso 


0147+400 


E 






100 


>26.6 


631 


563 


95 


~0.40 


scj?* 


0703+468 


E 






80 


>26.6 


1392 


1356 


75 


~0.32 


CSO 


0800+472 * 


E 






1000 


>26.8 


769 


447 


160 


~0.67 


scj?* 


0809+404 * 


G 




0.551 


1200 


26.90 


929 


522 


80 


0.28 


? 


0822+394 


G 




1.180 


50 


27.33 


1003 


942 


70 


0.30 


CSO 


0840+424A 


E 






80 


>26.8 


1243 


1187 


135 


~0.57 


cso 


1007+422 * 


E 






130 


>26.4 


370 


313 


265 


~1.12 


MSO 


1008+423 


E 






50 


>26.5 


521 


498 


115 


~0.48 


CSO 


1016+443 


G 


19.7 


0.33 R 


110 


26.11 


287 


266 


155 


0.44 


CSO 


1044+454 * 


G 


24.8 


4.10 K 


1000 


28.82 


352 


311 


175 


0.55 


? 


1049+384 * 


G 


20.9 


1.018 


100 


27.14 


574 


502 


210 


0.89 


CSO? 


1133+432 


E 






70 


>26.3 


1247 


1265 


45 


~0.20 


CSO 


1136+383 * 


E 






50 


>26.4 


400 


377 


70 


~0.28 


CSO 


1136+420 


G 


21.7 


0.829 


1000 


a 26.9 


402 


279 


120 


0.49 


? 


1159+395 


G 


23.4 


2.370 


50 


27.57 


535 


498 


70 


0.26 


CSO 


1225+442 


G 


18.2 


0.22 R 


200 


27.50 


316 


291 


420 


0.95 


CSO 


1242+410 


Q 


19.7 


0.811 


40 


27.08 


1215 


1126 


70 


0.29 


CSO? 


1314+453A* 


G 


21.8 


1.544 


160 


27.77 


574 


493 


185 


0.79 


CSO 


1340+439 


E 






70 


>26.5 


447 


397 


130 


~0.56 


CSO 


1343+386 * 


Q 


17.5 


1.844 


110 


27.70 


793 


717 


130 


0.54 


CSO 


1432+428B* 


E 






40 


>26.3 


809 


706 


50 


~0.21 


CSO 


1441+409 


E 






100 


>26.7 


834 


768 


130 


-0.56 


CSO 


1449+421 


E 






80 


>26.9 


639 


592 


120 


~0.50 


CSO 


2304+377 


G 


19.7 


0.40 R 


100 


26.44 


1299 


1217 


150 


0.49 


CSO 


2330+402 


E 






70 


>26.7 


730 


705 


100 


~0.43 


CSO 


2348+450 * 


G 




0.978 


200 


27.36 


630 


520 


310 


1.32 


MSO 


2358+406 


E 






80 


>26.8 


1159 


1161 


120 


~0.47 


CSO 



a the value given in Fanti et al. 2001 is wrong 

* steep-core Jet: elongated structure, very asymmetric in brightness, where the brightest component is at one source extremity 
and has a steep spectrum. 



additional presence of some weaker emission (bringing the 
overall size up to ~ 1", Table [j]) whose connection with 
the compact feature is unclear. In addition about half of 
the sample sources present a low frequency turnover in the 
total spectrum which is usually considered an indication 
of compactness in some sub-component. All these sources, 
therefore, deserved further investigation. 

The selected sample is presented in Table [TJ. The con- 
tent of the table is the following. 

Column 1 - Source name; an "*" refers to a note (Sect. |J); 

Column 2 - Optical Identification (Id) from Paper I (G = 
galaxy, Q = quasar, E = no known optical counterpart); 

Column 3 - R magnitude; 

Column 4 - redshift; "K" and "R" indicate that the red- 
shift is estimated by photometric measurements in the re- 
spective optical band (see Paper I for details); 

Column 5- VLA Largest Angular Size (LAS) in mas, from 
Paper I; 



Column 6 - log Po.4GHz (P hi W/Hz h~ 2 ); for E sources 
lower limits to the radio power have been computed as- 
suming z = 0.5 (see also Paper I for a wider discussion); 

Column 7 - flux density (mjy), interpolat ed at 1.67 GHz 
from the 1.4 GHz NVSS (Condon et al. |1998| ) and the 
5 GHz VLA data from Paper I; 

Column 8 - VLBA 1.67 GHz flux density (mjy), measured 
on the images presented here (see Sect. 3); 

Column 9 - VLBA Largest Angular Size (mas) ; 

Column 10 - Largest Linear Size (LLS, kpc h~ x ) from 
VLBA data; for E sources z = 1.05 was assumed (see 
Paper I for details) and LLS, given with a leading 
must be considered a lower limit; 

Column 11 - overall morphology (see Sect. |B|) 

The VLBA observations were carried out in February 
2000 at the frequency of 1.67 GHz, with a recording band- 
width of 32 MHz at 64Mbps. Each source was observed for 
a total of about 1 hour, spread into 10 to 15 short scans in 
order to improve the uv coverage. A typical example of uv 
coverage is given in fig. [TJ. As no single VLA antenna was 
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Fig. 1. Typical uv coverage 



available, the shortest baseline is Pie Town-Los Alamos 
(~1.3 MA). This may have caused some flux density losses 
in the more extended sources/components. 

3. Data reduction 

The data were correlated with the NRAO processor at 
Socorro. Amplitude calibration was obtained using the 
system temperature measured at the various antennas 
during the observations and the gain information pro- 
vided by NRAO. Consistency checks were performed with 
OQ208 and DA193, whose amplitudes were found to be 
generally within 2 % of what expected. 

Standard editing and fringe-fitting were applied. Most 
sources provided fringes with high signal-to-noise ratio on 
all baselines. 

Radio images were obtained using AIPS, after a num- 
ber of phase self-calibrations, ended by a final amplitude 
self-calibration. The last step was made with great care, 
as the gain self-calibration tends to depress the total flux 
density when extended components arc poorly sampled 
in uv coverage. Generally the gain corrections from self- 
calibration were < 3 %. In those sources in which they 
were > 3 %, or in which the total flux density decreased 
by more than 1 — 2%, they were not considered. 

The r.m.s. noise level (1 a), measured on the images 
far from the sources, is in the range of 0.1-0.2 mJy/beam, 
comparable to the expected thermal noise. The dynamic 
range, defined as the ratio of peak brightness to 1 a, is 
~ 700 on average, ranging from 300 to a few thousands. 
The typical resolution is ~ 4 x 8 mas. We must remark 
that the LAS determined on the VLBA images (Table 1, 
Col. 9) does not take into account the eventual presence of 



components completely resolved out by the VLBA obser- 
vations. We note also that it may exceed the VLA rough 
estimate due to the source structure (e.g. a double edge- 
brightened source will appear smaller when fitted by a 
single Gaussian with a resolution smaller than the sepa- 
ration between components). 

The total flux density "seen" by VLBA (Table 1, 
Col. 8) has been measured by integration over the source 
images using the AIPS task TVSTAT. The brightest sub- 
structures of each source were fitted by a gaussian model 
(task JMFIT) which allows to obtain total flux density, the 
beam-deconvolved Half Maximum Widths (HMW) and 
major axis position angle (p. a.). Sub-components are re- 
ferred to as North (N), South (S), East (E), West (W) 
and Central (Ce). When a component is split into more 
pieces, a digit (1,2, etc.) is added (e.g. N1,N2). Other 
symbols used are self-explanatory. The source images are 
shown in Fig. |^. The labels are marked on the images, next 
to each sub-component. The parameters of the brightest 
sub-components are given in Table which has the fol- 
lowing layout: 

Columns 1,2- Source name and sub-component label; 
an "*" indicates that the component is extended and its 
flux density is obtained by TVSTAT (see text); 

Column 3 - VLBA flux density; 

Column 4 - deconvolved angular sizes of major and minor 
axis of a Gaussian component and position angle of the 
major axis as estimated using the AIPS Gaussian-fitting 
task JMFIT. 

For more extended features, marked by an "*" in 
Table 2, flux densities are evaluated by means of TVSTAT. 
For extended structures underlying more compact ones, 
the flux density is obtained as difference between the in- 
tegrated flux density and the sum of the sub-component 
gaussian fit flux densities. 

We have checked the integrated flux densities in our 
images with the low resolution flux densities obtained by 
interpolating to our observing frequency the values from 
the NVSS (Condon et al. |1998[ ) and the from VLA data 
at 5 GHz from Paper I. The comparison is displayed in 
Fig. | 

There is clearly a systematic flux density difference, 
typically J> 5 %, which could be partially ascribed to cali- 
bration uncertainties. For the more extended objects there 
are larger flux density losses, clearly due to the absence of 
short uv-spacings (see, as an example, Fig. 0). The two 
sources with the largest deviations in Fig. g (0800+472 
and 0809+404) do possess indeed another component de- 
tected by the VLA but completely resolved out by the 
VLBA baselines (see also next section). 

4. Comments on individual sources 

- 0800+472: the source has an overall extent of ~ 1 arcsec 
in the VLA 8.4 GHz image (see paper I). About 40 % 
of the flux density is missing in our VLBA image. Since 
~ 90 % of the 8.4 GHz flux density originates within a 
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total flux density (Jy) 



LAS (mas) 



Fig. 2. (left): Comparison between the VLBA integrated flux densities and total flux densities interpolated from the 
literature. The line has a slope of one. (right): the ratio between the VLBA and the total flux density is shown as a 
function of the overall source angular size 



region of ~ 120 mas, it is likely that most of the missing 
flux density is on angular scales of this order. The VLBA 
source structure is strongly bent and the bending probably 
continues to match the arcsec structure. It is possible that 
we are in presence of an asymmetric core-jet, with a weak 



core, as in Stanghellini et al. (2001) 



- 0809+404: at arcsec scale the source appears as a very 
asymmetric double, with 1.2" separation and a flux den- 
sity ratio ~ 100:1 at 4.9 GHz (see paper I). The structure 
we see here corresponds to the brighter component and 
is rather amorphous, with the p. a. of the major axis mis- 
aligned by ~ 40° with respect to the orientation of the 
outer weak component. About 45 % of the total flux den- 
sity is missing in our VLBA observations. 

- 1007+422: about 15 % of total flux density is missing in 
our VLBA observations. 

- 1044+454: on the arcsec scale the source appears as a 
very asymmetric double, with flux density ratio ~ 20:1 at 
4.9 GHz (see paper I). On the VLBA scale the structure 
of the northernmost brighter component is again double 
and asymmetric, along the same p. a. but on a ~ six times 
smaller scale. The source could be either a core-jet with a 
weak core or a very asymmetric MSO. About 12 % of the 
total flux density is missing in our VLBA observations 

- 1049+384: About 12 % of the total flux density is missing 
in our VLBA observations 

- 1136+420: on the arcsec scale the source appears as a 
very asymmetric double, with flux density ratio ~ 20:1 at 
4.9 GHz (see paper I), roughly oriented in the East-West 
direction. At VLBA resolution only the brighter compo- 



nent is visible and displays a very amorphous structure. 
About 30 % of the total flux density is missing in our 
VLBA observations 

- 1314+453: about 15 % of the total flux density is missing 
in the VLBA image. 

- 1343+386: a weak feature, accounting for ~ 9 mJy is 
present to the North of the major components. 

- 1432+428: a weak feature, accounting for ~ 10 mJy 
is present to the South of the major components. About 
13 % of the total flux density is missing in our VLBA 
observations 

- 2348+450: about 20 % of total flux density is missing in 
the VLBA image. 

5. Discussion 

The VLBA observations have produced images with suffi- 
cient resolution to describe the source morphology for the 
large majority of cases. 

The present observations do not allow to pick up the 
radio core, due to the relatively low observing frequency, 
but the overall morphology strongly indicates that the 
sources are mostly two-sided. We see either well sepa- 
rated lobes or lobes connected by bridges, possibly jets. 
Most sources are edge brightened, suggesting the pres- 
ence of hot spots, in general not fully resolved from the 
mini-lobe emission, which contribute most of the source 
flux density. We tentatively classify these sources as CSOs 
(two of them are actually MSOs). There are also a few ob- 
jects with edge darkening, as B3 2304+377, B3 2348+450 
and B3 2358+406. There are a number of sources where a 




Fig. 3. VLBA images: the first contour (f.c) is generally three times the r.m.s. noise level on the image; contour levels 
increase by a factor of 2; the restoring beam is shown in the bottom left corner in each image. 
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Fig. 3. VLBA Images (cont.) 
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Fig. 3. VLBA Images (cont.) 
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Fig. 3. VLBA Images (cont.) 
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Fig. 3. VLBA Images (cont.) 



compact component, centrally located, could be the source 
core. However, having only one relatively low frequency 
we are not in the position to make a definitive statement 
about core identifications. 

Some asymmetries in the flux densities of the two sides 
are seen, as also some moderate distortions. B3 0800+472 
and B3 2304+377 are extreme examples. The sources 
B3 1225+442 and B3 1449+421 are reminiscent of an 5 



morphology (see also Stanghellini et al. 2001 ). 

The average ratio between the component size trans- 
verse to the source axis (#trans) and Largest Angular Size 
(LAS) is 9 tIsns /LAS ~ 0.08 ± 0.02, consistent with find- 



ings in other source samples (e.g. Fanti et al. 1990) 

We have computed the equipartition parameters for 
the source components, under the following assumptions: 



a) proton to electron energy ratio of one; b) filling factor 
of one; c) maximum and minimum electron (and proton) 
energies corresponding to synchrotron emission in the fre- 
quency range 100 GHz to 10 MHz; d) ellipsoidal volumes 
with axis corresponding to the observed ones. We also 
computed the brightness temperatures (Tq) of each com- 
ponent. In order to make easier the reading of the paper, 
we do not report the individual component values, since 
they can be obtained from the parameters in Table but 
only mention the typical values. 

The magnetic fields are in the range of several mGauss, 
going up to ~ 20 mGauss in the most compact components 
(e.g. B3 1133+432). Brightness temperatures range from 
- 10 8 up to < 10 11 K. 
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Table 2. Observational parameters of source components 



Source 




s 


Si 


2 


(u a) 


Source 




s 


Si 


g 2 


(u a) 






(m Jy) 












(ni Jy) 








I'D 


(2) 


(3) 




(4) 






(2) 


(3) 




(4) 


_ 


0039+373 


Nl 


229 


4.7 


4.3 


101 


1225+442 


El 


138 


26.6 


15.1 


88 




N'2 


143 


11.0 


6.8 


17 




E* 


28 










N* 


17 










Wl 


94 


46.6 


23.3 


71 




SI 


303 


13.2 


6.9 


24 




W2 


35 


32.4 


10.6 


134 




s* 


88 








1242+410 


S 


818 


20.8 


5.1 


26 


0147+400 


E 


208 


2.7 


2.4 


41 




N 


208 


7.7 


4.8 


53 




Wl 


134 


11.6 


7.6 


40 




T* 


101 










W2 


109 


13.7 


10.6 


26 


1314+453 


Wl 


175 


26.7 


15.6 


71 




T* 


113 










W2 


77 


31.6 


26.9 


73 


0703+468 


W 1 


794 


7.1 


4.01 


51 




Ce 


60 


19.0 


15.0 


143 




w* 


22 










El 
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28.0 


14.6 


63 




El 


497 
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4.0 
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E2 


76 
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10.2 
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E* 


44 










NW 


18 
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82 


0800+472 


SI 
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49 
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60 
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24 
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11.3 
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S 


27 


12.0 


6.6 


83 


0809+404 


N 
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16.1 


12.2 


161 


1343+386 


SI 
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5 . 9 
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161 
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6. Conclusions and future plans 

We have presented the results of VLBA observations of 
the most compact sources (LAS <^ 1") in the B3-VLBA 
CSS sample from paper I. 

The large majority of sources have a "double type" 
structure, with sizes ranging from ~ 200 h^ 1 pc to ~ 1 h^ 1 
kpc. Together with the VLA data on the large size sources 
and wi th EV N & MERLIN data, presented in Dallacasa 
et al. (|2002| , Paper HI) on intermediate size sources, we 



have at the end a new determination of the Linear Size 
distribution in the range ~ 0.2 hr 1 < LS'(kpc) < 20 hr 1 . 
This is an important piece of information to test the source 
evolution model. We will discuss it in a separate paper. 

The sources presented here have just been observed 
with the VLBA at higher frequencies, with the aim of de- 
tecting the radio core, therefore obtaining an unambigu- 
ous classification, and of studying the structure of the hot 
spots. The results will be presented in a forthcoming pa- 
per. 
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Fig. 4. Plot of the correlated flux density versus projected 
baseline for the source B3 1314+453. The horizontal ar- 
rows indicate the integrated flux density in the VLBA im- 
age and total flux density from the VLA as given in Table 
1 
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